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1.  Introduction 

The  paper  on  shock  waves  In  elastic-plastic  solids  [1]*  by  Germain  and 
the  author  of  this  report.  Is  concerned  In  the  main  with  the  application  of 
a generalization  of  classical  rate-independent  constitutive  relations  for 
elastic-plastic  media  to  the  analysis  of  the  propagation  of  plane  compressive 
waves.  Waves  exhibiting  both  continuous  and  discontinuous  stress  profiles  were 
considered.  Since  the  generation  of  shock  waves  is  Intimately  associated  with 
non-linear  effects,  the  generalization  of  classical  theory  adopted  Incorporated 
finite  deformation  kinematics.  It  was  found  possible  to  formulate  the  theory 
so  that  the  mean  hydrostatic  pressure  component  of  the  stress  tensor  and  the 
deviator  or  dlstortlonal  components  appeared  In  separate  terms  as  did  the  dila- 
tation and  the  strain  deviator.  Thus  the  extensive  literature  on  pressure- 
volume  relations  developed  In  connection  with  the  hydrodynamic  theory  of  shock 
waves  In  metals  could  be  utilized  directly  for  the  analysis  of  wave  propagation 
In  metalswhen  shear  strength  also  influences  the  motion.  The  separation  of 
pressure  and  shear  effects  was  carried  out  on  the  basis  of  the  structure  of  the 
free  energy  function  for  thermo-elastlc-plastlc  deformation  In  a physically 
meaningful  manner  which  also  satisfied  the  invariance  requirements  of  nonlinear 
continuum  mechanics.  For  example,  the  latter  links  the  influences  of  density 
change  on  direct  stress  and  on  shear  stress  components  which  follows  from  the 
structure  of  finite  deformation  thermo-elastic  theory. 

In  the  analysis  of  shock  waves  Involving  plastic  flow  using  rate  Independent 

theory  [1],  it  was  found  that  a shock  wave  solution  could  not  be  determined  on 

this  basis  when  work  hardening  was  considered  or  %rhen  part  of  the  plastic  work 

expended  la  transformed  into  Internal  energy  associated  with  breakdom  of  the 
_ 

Numbers  In  square  brackets  refer  to  the  references  listed  at  the  end  of  the 

report. 
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crystal  lattice  by  the  generation  of  dislocations.  This  component  of  energy 
absorption  by  the  metal  has  a negligible  Influence  on  the  thermo- elastic  state 
of  the  element  and  so  does  not  Increase  the  entropy.  In  contrast  to  the  rest 
of  the  plastic  work  which  Is  dissipated  Into  heat.  Thus  knowledge  of  the 
work  done  on  the  material  during  passage  of  a shock  wave  does  not  determine 
the  new  Internal  energy  In  the  thermo-elastic  system  since  part  Is  absorbed 
In  forming  the  dislocation  distribution.  This  prevents  direct  evaluation  of 
the  state  after  passage  of  the  wave.  Similarly  the  work  absorbed  In  plastic 
flow  which  determines  the  yield  stress  of  the  work-hardening  metal  after  passage 
of  the  shock  wave  can  only  be  calculated  from  a study  of  the  shock  structure 
which  calls  for  a rate  type  plasticity  law  to  correctly  Incorporate  the  energy 
dissipated  In  the  shock  wave.  Plastic  work  Is  expressed  by  a time  Integral  of 
a product  of  the  yield  stress  and  the  plastic  strain  rate  and  cannot  be  deter- 
mined directly  In  terms  of  variables  only  expressing  properties  before  and  after 
passage  of  the  wave. 

In  view  of  these  difficulties,  a hypothetical  rate-dependent  plasticity 
law  was  postulated  which  permitted  the  development  of  a theory  which  would 
determine  the  shock  structure  for  a steady  wave  and  hence  the  variations  of 
strain,  temperature,  plastic  strain  and  plastic  work  with  position  through  the 
wave  and  hence  the  thermo-elastic  state  of  the  material  after  passage  of  the 
wave.  Moreover  it  was  shown  that  the  choice  of  rate-law  for  the  plastic  strain- 
rate  did  not  influence  the  sequence  of  values  taken  on  by  these  dependent  vari- 
ables but  only  their  profiles  In  space.  The  plots  of  any  three  of  the  dependent 
variables  e.g.  temperature,  p^  itlc  strain  and  plastic  work  against  the  other, 
i.e.  the  strain,  would  be  Independent  of  the  rate  law.  In  spite  of  this.  It  Is 
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worthwhile  to  consider  experimentally  motivated  rate-laws  so  that  actual  wave 
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profiles  can  be  evaluated.  The  results  mentioned  concerning  the  deduction  of 
partial  Information  about  the  solution,  without  knowledge  of  the  precise  rate 
law,  Is  a generalization  of  the  Rayleigh  line  concept  for  steady  waves  In  ma- 
terials defined  by  simpler  constitutive  relations.  For  example.  In  the  theory 
of  shock  waves  In  gases,  or  equivalently  hydrodynamic  analysis,  the  knowledge 
of  the  work  done  on  the  material  due  to  passage  of  a shock  Is  sufficient  to 
determine  the  final  thermodynamic  state  of  the  material.  Essentially  because 
of  the  additional  Internal  variable:  plastic  strain,  the  more  complicated  situ- 
ation already  described  arises  in  plastic  analysis. 

In  the  next  section  the  formulation  of  thermo-elastlc-plastlc  constitutive 
relations  for  metals  Incorporating  a rate  of  plastic  strain  Influence  will  be 
discussed,  and  compared  with  the  rate- Independent  form.  Subsequently  the  appli- 
cation of  these  for  analyzing  plastic  wave  propagation  phenomena  will  be  con- 
sidered. In  discussion  of  the  Influence  of  rate-dependent  analysis,  other 
characteristics  which  generate  analogous  effects  on  wave  profiles  will  be  con- 
sidered. ^ 

2.  Finite  Deformation  Thermo-Elastlc-Plastlc  Theory 

The  constitutive  relations  to  be  applied  comprise  a generalization  of  the 
finite  deformation  theory  described  In  reference  [2].  Because  non-linear  effects 

t 

exert  a major  Influence  on  shock  wave  generation  and  propagation,  It  Is  Impor- 
tant to  have  a correctly  formulated  constitutive  relation  which  Incorporates 
both  geometric  and  material  nonlinearities.  The  kinematics  [2]  Is  based  on  the 
matrix  product  representation  of  the  combined  effect  of  elastic  and  plastic 
deformation: 

F - f'e”*  , (1) 


where  Is  the  deformation  gradient  matrix  3x^/3Xj  where  x^  is  the 

position  vector  of  an  element  in  the  deformed  configuration  and  for 

the  element  in  its  reference  configuration.  is  the  deformation  gradient 

matrix  for  plastic  deformation  relative  to  the  initial  undeformed  state  and 
F^  the  elastic  deformation  gradient  relative  to  a state  unstressed  after 
plastic  deformation.  At  the  current  time  t,  F - F (X«t)  since  in  general 
the  deformation  of  a body  will  be  distributed  non-homogeneously  . 

For  the  study  of  the  propagation  of  plane  waves  of  one-dimensional  strain 
as  generated  in  a plate  slap  experiment,  principal  directions  remain  fixed  in 
space  and  in  the  body  normal  to  and  in  the  plane  of  the  wave  front.  Thus  the 
principal  components  only  of  F need  to  be  Introduced: 


Ih  0 0 

F - ( 0 X2  0 

\o  0 X. 


Components  X^  are  the  stretch  ratios  of  the  length  In  the  direction  i after 
and  before  the  deformation,  and  X^  = X^  for  isotropic  media  because  of  sym- 
metry in  the  plane  of  the  wave.  Then  taking  natural  logarithms  of  the  total, 
elastic  and  plastic  stretch  ratios  determines  the  natural  strains: 


- In(X^)  , e*  - ln(X*)  , - In  (X^) 


and  the  relation  (1)  for  total  deformation  reduces  to  additivity  of  the  elastic 
and  plastic  strain  components 


e p 

+ Cj 
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as  discussed  In  [3]  and  [1].  For  plane  wave  analysis  the  strains  will  In 
general  be  functions  of  time  t and  of  position  along  an  axis  normal  to  the 
wave-front.  The  latter  could  be  the  Lagrange  coordinate  In  the  reference 

geometry  or  the  Euler  coordinate  In  the  current  configuration. 

The  stress  o and  entropy  s are  given  In  terms  of  the  density  p , the 
the  elastic  deformation  and  temperature  9 through  the  Helmholtz  free 

energy  function  ¥ by  the  thermo-elastic  relations: 


p F®3¥/3F® 

(5) 

- 3¥/3e 

(6) 

As  discussed  In  [4]  and  [5]  the  free  energy  function  Is  taken  to  have  the 
form 


¥ - ¥ (F®,0)  +¥(w) 


(7) 


where  w Is  the  plastic  work  per  unit  mass.  The  first  term  on  the  right 
hand  side  expresses  the  thermo-elastic  response  and  the  second  term  the 
energy  stored  In  the  dislocation  distribution  resulting  from  plastic  flow. 
The  latter  has  little  Influence  on  thermo-elastlclty  as  confirmed  by  the 
Insensitivity  of  elastic  constants  to  plastic  flow.  The  free  energy  com- 
ponent ¥ constitutes  only  a small  part  of  the  plastic  work  w 


¥ - (1-y)  w , Y “ 0. 9 


(8) 


the  rest  being  dissipated  Into  heat. 

The  thermo-elastic  component  of  the  free  energy  Is  taken  to  have  the 

form 
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(9) 


't  (e®  , E®  , E®  , e)  « f (Ij^.e)  + h (J2,J3) 


where  I.  is  the  first  invariant  of  the  elastic  strain  tensor: 


e , e 
ei  + E2 


+ E, 


(10) 


and  J2  and  are  the  second  and  third  invariants  of  the  elastic  strain 

devlator  tensor  e 


■^2  " 2 

J3  “ 3 [(e®)^+  (ep^  + (e®)^l 


(11) 

(12) 


Since  is  equal  to  the  volumetric  strain,  plastic  flow  being  incom- 

pressible, f(I^,6)  prescribes  the  pressure,  volume,  temperature  relation 
of  the  metal.  Under  general  loading  the  pressure  p is  defined  as 
- « the  average  of  the  normal  compressive  stress  components,  and  the 

same  dilatation  law  arises.  The  function  h(J2«J3)  prescribes  the  deviator 
or  distortional  components  of  the  stress  and  has  been  further  simplified  to 
2iiJ2/PQ  In  the  analysis.  The  temperature  6 is  not  included  in  this  term 
since  shear  strain  does  not  usually  stimulate  thermo- elastic  coupling. 

In  plane  wave  experiments,  prior  to  the  arrival  of  release  waves  from 
the  free  lateral  surface  of  the  specimen,  the  lateral  strains  are  zero: 


C2  - £3  - 0 (13) 

and  since  plastic  deformation  is  incompressible 

C3  - - In  (Pq/p)  (14) 

where  p^  is  the  density  In  the  reference  state. 
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For  this  special  deformation  the  stress  components  can  be  expressed  in  terms 


of  the  longitudinal  total  and  plastic  strain  components,  and  respec- 
tively, by  using  (4)  and  (13)  to  obtain  elastic  strains  and  then  applying  (5) 
[see  [1],  eq.(44)] 


o - -p+iiL£e  - 

0 -a  .-p-iPPe  +iie.  -P 
2 3 3 Pq  1 Pq  1 


The  pressure  p is  given  by  the  hydrodynamic  relation  already  mentioned  and 
prescribed  by  the  function  f(lj^,0).  It  is  often  convenient  to  express  the 
pressure  - temperature  - dilatation  relation  in  the  form 


PCs.Ej^)  , 0(s,c^) 


since  s can  be  a convenient  Independent  variable  for  adiabatic  loading 

Equations (15)  and  (16)  indicated  that  the  separation  of  dilatation  and 
deviator  terms  in  (9)  does  produce  a coupling  of  the  deviator  stress  com- 
ponents with  volume  change  according  tothe  factor  P/Pq  • 

We  note  that  in  addition  to  the  total  longitudinal  strain  c^,  the 
plastic  longitudinal  strain  component  c^  acts  as  an  additional  state  vari- 
able which  permits  determination  of  the  stress  for  prescribed  temperature  or 
entropy.  This  is  so  without  any  statement  concerning  the  laws  of  plasticity. 

The  rate  Independent  theory  of  plasticity  described  in  detail  in  [2]  and 
[1]  adjoins  a work-hardening  yield  condition  which  must  be  satisfied  for 
plastic  flow  to  take  place,  in  the  form  that  an  isotropic  function  of  the 
deviator  of  p ^o  is  prescribed  by  the  plastic  work  w and  the  temperature 


0 . For  our  particular  problem  this  condition  is  contained  in  the  expression: 
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Y(w,0)  + PgP~^(Oj^-  o^)  1 0 


where  Y is  the  yield  stress  in  tension  or  compression  at  the  initial 
density  and  temperature  6 ; the  upper  sign  corresponds  to  longitudinal 
tensile  yield  and  the  lower  to  compressive.  The  equality  sign  permits 
plastic  flow  to  take  place  and  the  inequality  indicates  that  the  element  is 
inside  the  yield  surface  and  hence  currently  subject  only  to  elastic 
strain  Increments.  When  plastic  flow  takes  place  the  plastic  work  per 
unit  mass  is  determined  from: 


or,  using  (18) 


Pq  w “ ± Y (w,0) 


with  the  same  sign  convention. 

The  relations  prescribed  are  appropriate,  for  example,  to  determine 
the  stresses  generated  by  a continuous  adiabatic  process  defined  by  a 
given  variation  of  e^(t)  from  a prescribed  state.  If  the  element  is  ini- 
tially subjected  to  stress  inside  the  yield  surface,  so  that  the  Inequality 
in  (18)  is  satisfied,  only  elastic  Increments  of  strain  can  occur  so  that 


.p  • • 

Cj  • w - s “ 0 


and  (15)  and  (16)  with  p determined  by  (17)  give  the  stress  variation 
since  is  known  for  the  initially  prescribed  state.  When  plastic  flow 

takes  place  the  equality  (18)  %rith  (15)  and  (16)  yield 


I 


(22) 


3eP/2  - + Y/(2p) 


The  part  yv  of  the  plastic  power  which  Is  dissipated  Into  heat  causes 
an  entropy  Increase  according  to 


9 s “ Y w 


(23) 


Combining  this  with  (20)  and  changing  the  Independent  variable  from 
yields  the  differential  equation  system 


Q ds 

D 0 — 

^0  de. 


dw 


t to 


(24) 


The  variable  e^  can  be  eliminated  using  (22),  and  utilizing  the  known 
functions  0(s,ej^)  from  (17)  and  Y(w,0)  gives  a pair  of  differential 
equations  to  determine  w(c^)  and  s(c^)  using  the  known  Initial  state 
to  provide  Initial  conditions.  The  plastic  strain  component  Is  given 

by  (22),  Y(w, 0)  being  obtained  using  0(s,c^),  and  then  the  stresses  can 
be  evaluated  from  (15)  and  (16).  Note  that  the  time  variable  does  not 
appear  explicitly  In  the  solution  for  either  elastic  or  elastic-plastic 
deformation,  in  conformity  with  the  adoption  of  the  rate  Independent  plas- 
ticity constitutive  relation.  It  might,  for  example,  be  convenient  to 
utilize  the  variable  as  a pseudo  "time"  In  all  rate  relations  considered 

so  far  In  this  section,  for  example  equations  (19)  and  (20),  and  evaluate  the 
stresses  without  Introducing  real  time  t . 

Consider  now  a modification  of  the  hypothetical  rate  dependent  plasti- 
city law  of  [1]  to  bring  It  Into  conformity  with  experimental  findings. 
Equations  (15)  and  (16)  show  that  the  stress  can  be  determined  In  terms  of 


total  and  plastic  longitudinal  strain  components.  These  are  deduced  from 
the  thermo-elastic  stress  law  (5)  after  expressing  the  elastic  strain  com- 
ponents In  terms  of  and  by  means  of  (4),  (13)  and  the  incompressi- 
bility of  plastic  flow.  These  stresses  can  be  substituted  Into  the  yield 
condition  and  a rate  of  plastic  strain  prescribed  In  terms  of  the  extent 
to  which  the  yield  condition  has  been  violated.  If  the  stress  point  lies 
Inside  or  on  the  yield  surface  the  rate  of  growth  of  plastic  strain  will 
be  zero. 

Whatever  the  law  of  plasticity  (15)  and  (16)  yield  the  relation: 


3 -0 

^1  2 1 “ p 2ii 


so  that  the  yield  condition  (18)  can  be  written: 


(25) 


Y(w,e)  + 2p(e^-  I eP)>  0 


(26) 


Thus  a convenient  means  of  expressing  over-stress  Is  through  the  variable 
z , defined  by 


/ 3 p,2 

’a? 


(27) 


This  will  be  positive  when  the  yield  condition  Is  violated  and  negative 
If  the  stress  point  lies  inside  the  yield  surface.  Thus  an  over-stress 
condition  of  the  Malvern  type  can  take  the  form 


(^1  - f E^)  k(z) 


(28) 


where  k Is  an  appropriate  function.  The  factor  (e^  - e^)  is  Intro- 
duced to  ensure  the  correct  sign  of  the  deduced  plastic  strain  rate.  The 
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yield  stress  function  Y(w,6)  now  has  a changed  significance  since  at  yield 
the  plastic  strain-rate  will  be  zero,  and  the  yield  stress  must  be  exceeded 
for  non-zero  plastic  strain-rate  to  occur.  Thus  k is  a monotonlcally 
increasing  function  with  k(0)  ■ 0. 

Using  (25)  the  rate  of  plastic  work  is  now  given  by 


w 


pQ  ^*^1  2 ^1^  ^1 


(29) 


The  adiabatic  entropy  growth  relation  (23)  can  still  be  expected  to  apply,  and 
(28),  (29)  and  (23)  combined  with  the  thermo-elastic  relation  specify  the  corres- 
ponding thermo-elastlc-plastlc  constitutive  relation. 

If  a prescribed  variation  e^(t)  from  a known  state  is  imposed  under  adia- 
batic conditions, (23) , (28)  and  (29)  constitute  a system  of  three  differential 
equations  to  determine  w , and  s , when  use  is  made  of  the  known  function 
6(s,ej^)  from  (17). 

The  introduction  of  a strain-rate  effect  in  combination  with  a work  har- 
dening law  Introduces  "a  strain-rate  history  effect"  of  the  type  demonstrated 
by  Klepaczko  [6]  and  mentioned  by  Clifton  [7]  (p.l03)  in  his  article  on  plastic 
waves.  Thus,  for  a given  plastic  strain  Increment,  more  plastic  work  is  expended 
if  the  straining  is  carried  out  at  higher  strain-rate  since  the  stress  is  higher 
according  to  (28),  (27)  and  (25).  If  straining  is  first  carried  out  with  an 
initial  strain-rate,  , until  a prescribed  strain  is  achieved,  and  thereafter 

s 

the  strain-rate  is  changed  to  ; for  a particular  strain  during  the  second 
part  of  the  stepped  strain- rate  test,  the  stress  would  be  higher  had  the  initial 
strain-rate  been  bigger,  since  the  plastic  work  w and  hence  the  equili- 

brium yield  stress  Y(w,6)  would  be  larger.  Qualitatively,  such  behavior  has 


BP 


been  observed  by  Klepaczko  [6].  As  discussed  In  the  following  section,  such  a 
strain-rate  history  dependent  type  law  Is  Incorporated  Into  the  analysis  of 
shock  wave  structure  presented  In  [1]. 

By  taking  the  yield  stress  Y to  depend  on  , Instead  on  w , 
which  In  the  case  considered  In  this  report  of  longitudinal  monotonlc 
straining  Is  equivalent  to  the  commonly  used  generalized  plastic  strain 
? . dt)  (see  Hill  [8]  p.  30),  the  appearance  of  a rate 

history  effect  Is  avoided,  since  then  the  rate  at  which  the  stress  history  Is 
carried  out  does  not  affect  the  equilibrium  yield  stress  Y . As  described 
In  the  following  section  this  change  does  not  greatly  modify  the  analysis  of 
straining  In  wave  propagation  based  exactly  on  the  elastic-plastic  theory 
presented,  although  a good  approximation  Is  more  easily  obtained.  A some- 
what similar  plastic  strain-rate  and  total  plastic-strain  dependent  thermo- 
elastlc-plastlc  constitutive  relation  has  been  discussed  by  Clifton  [7](p.l31 
ff.).  This  also  contained  a vector  of  parameters  which  characterized  the 
Internal  structure  of  the  material.  The  plastic  work  w could  be  such  an 
Internal  variable  and  the  work-hardening  theory  can  be  considered  to  be 
analogous  to  such  a law. 

3.  Wave  Propagation  Phenomena 

As  mentioned  In  the  Introduction,  certain  difficulties  In  obtaining  shock 
wave  solutions  for  strain-rate  Independent  elastic-plastic  materials  were 
encountered  for  work-hardening  plasticity  laws.  Related  rate-dependent  laws 
were  therefore  studied  which  result  In  dual  shock  structure  of  an  elastic  wave 
exhibiting  discontinuous  properties  attached  to  a steady  wave  associated  with 
tile  development  of  plastic  flow.  This  corresponds  to  an  over-driven  wave  forced 
to  travel  at  a speed  higher  than  elastic  wave  velocity  (see  [7],  p.  148  or 
Herrmann  [9],  p.  10).  It  la  not  uncommon  to  refer  to  the  combination  as  a 
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shock-wave,  terming  the  variations  of  properties  in  the  steady  wave  the 
study  of  shock  structure.  We  shall  term  these  two  components  the  shock  wave 
and  the  steady  wave  with  invariant  profile. 

A number  of  peripheral  questions  arise  in  the  theoretical  analysis  of 
the  phenomenon  and  in  connection  with  experimental  measurements.  For  example 
it  was  commonly  the  practice  to  analyze  steady  waves  as  having  an  invariant 
profile  for  measurements  taken  at  a fixed  point  in  space,  that  is  to  measure, 
say,  the  velocity  v(xj^,t)  for  varying  t and  fixed  x - x^^  . The  Euler 
coordinate  x^  denotes  a fixed  point  in  space  as  the  material  moves  across 
it.  It  is  now  realized  that  measuring  devices  are  commonly  fixed  to  a material 
element  and  hence  yield  the  data  v(X^,t)  where  the  Lagrange  coordinate 
defines  a position  in  the  reference  configuration,  that  is  at  a particular 
material  element.  It  seems  plausible  that  a wave  which  has  an  invariant 
profile  at  each  point  in  space  would  also  have  one  at  each  point  in  the  body. 
However,  since  the  constitutive  relation  for  elastic  plastic  media  is  highly 
non-linear,  it  seems  unlikely,  without  further  consideration,  that  the  trans- 
formation from  one  coordinate  to  the  other  could  maintain  the  invariance. 

The  motion  along  the  axis  x of  the  cross-section  labelled  X in  the 
reference  configuration  in  Fig.  1 is  given  by 


X ■ x(X,t) 


(30) 


where  the  subscript  to  the  coordinates  has  been  dropped  for  one-dimensional 
waves.  The  stretch  ratio  X is  given  by 
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Hence 


. 3x 
^ “ 3X 


li  a X 3v 
3t  * 3t3X  “ 3X 


where  v is  the  particle  velocity 


If  there  is  a steady  wave  in  the  Lagrange  frame  which  travels  without 
change  of  profile  when  considered  plotted  in  the  reference  configuration, 
then  dependent  variables  a,v  and  X are  functions  of  (X  - only, 

when  is  a constant  - the  wave  speed  in  the  Lagrange  sense,  l.e.  along 

the  reference  configuration.  Thus 


X - X(X-Cj^t)  , V = v(X-Cj^t) 


Relation  (32)  then  takes  the  form 


- C.  X'  ■ v' 

w 


where  the  prime  denotes  differentiation  with  respect  to  the  single  argument. 
Integrating  (35)  gives 


1 


V + C, X • constant 


(36) 


Now,  corresponding  to  propagation  velocity  C along  the  reference  frame, 

L 

the  propagation  velocity  along  the  deformed  body  (l.e.  in  space)  is 

Cj^X  + V (37) 

since  travel  dX  in  the  reference  state  is  travel  relative  to  the  deformed 
material  through  the  distance  AdX  , the  material  meanwhile  moving  with 
velocity  V . This  is  the  Euler  wave  speed 

Cg  - Cj^X  + V (38) 

and  by  (36)  is  therefore  constant.  Since  constant  values  of  X,v  and  a 
propagate  with  velocity  C in  the  reference  frame  tilth  coordinate  X , 

la 

they  travel  with  Cg  in  space,  so  that 

X ■ X(x  - Cgt) , - etc.  (39) 

Thus  (34)  corresponds  to  s steady  wave  in  space  (Euler  coordinates). 

Note  that  for  a steady  wave  traveling  into  an  undisturbed  body  (e.g. 

as  discussed  in  [7]  p.  148  and  [9]  p.  10)  X - 1 and  v > 0 at  the  wave 

front  so  that  C-  ■ C by  (38).  Of  course  the  profiles  of  stress  o , 

E L 

velocity  v or  stretch  X are  not  the  same  expressed  as  functions  of  X and  x 
Since,  for  a steady  wave  to  develop,  it  must  be  propagating  into  material 
homogeneously  stressed  and  at  uniform  velocity,  if  axes  are  chosen  which 
bring  the  material  ahead  of  the  wave  to  rest,  and  if  material  uniformly 
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strained  as  that  ahead  of  the  wave  is  used  as  the  reference  state,  the 
Euler  and  Lagrange  wave  speeds  will  be  the  same  for  all  cases. 

To  extend  the  consideration  slightly,  following  Clifton's  development 
[7]  p.  149),  Che  equation  of  motion  In  the  reference  frame  Is 


3a  ^ 3v 

3x  ‘‘oat 


(40) 


Since  there  is  no  lateral  motion  In  one-dimensional  strain  situations, 
cross-sectional  areas  are  Identical  In  the  reference  and  current  configurations 
so  that  the  nominal  stress  in  (40)  Is  equal  to  the  usual  Cauchy  or  "true" 
stress  In  the  deformed  body.  For  a steady  wave  with  velocity  C^^  In  the 
reference  frame;  (40)  becomes 


o'  • - Po  (41) 

which  combined  with  (35)  gives 

o’  - PqCJ  X'  (42) 

which  Integrated  Is: 

2 

a ■ p-  C^  X + constant  (43) 

0 L 

Thus  the  stress  and  stretch  ratio  are  linearly  related  In  a steady 
wave,  this  relationship  cosssonly  being  called  the  Rayleigh  line.  Thus  In 
spite  of  the  non-linear  Irreversible  constitutive  relation  governing  rate 
dependent  clastic-plastic  response  at  finite  strain,  Che  constraint  of  a 


steady  wave  forces  such  simple  connection  between  the  variables.  It  is 
Important  to  observe  that  these  results  depend  only  on  the  compatibility  of 
the  kinematic  variables  (32)  and  the  equation  of  motion  (40)  and  their  form 
does  not  depend  on  the  constitutive  relation  of  the  material. 

It  is  thus  clear  that  these  results  follow  directly  from  the  conserva- 
tion laws  as,  for  example,  developed  in  [1]  Section  2.  The  discussion  of  this 
aspect  is  related  to  the  analogous  studies  by  Duvall  [10]  and  Johnson  and 
Barker  [11].  For  the  x-axis  chosen  so  that  the  wave  is  propagating  in  the 
direction  of  increasing  x into  material  at  rest,  this  coordinate  being  de- 
noted by  X , let  v(x,t)  be  the  particle  velocity.  The  relation  for  conser- 
vation of  mass  for  a wave  of  one-dimensional  strain  is  ([1],  section  2): 

1^  + ^ (pv)  - 0 (44) 

If  the  Euler  wave  velocity,  C , reversed,  thus  now  directed  towards 

E 

decreasing  x , is  superposed  on  the  whole  configuration,  the  modified  x 
coordinate  will  be: 


X - X - Cgt  (45) 

and  the  particle  velocity  relative  to  this  coordinate  system  will  be: 

V ■ V - C (46) 

E 


Tba  steady  wave  will  be  transformed  into  a stationary  wave,  so  that 


The  other  dependent  variables  p,o  and  X will  not  change  but  will  become 
functions  of  x only  for  the  stationary  wave: 

p(x,t)  - p(x),  a(x,t)  - o(x),  X(x,t)  « X(x)  (48) 

The  equation  for  conservation  of  mass  will  now  be 

(pv)  “ 0 (49) 

which,  since  p > p(x),  reduces  to 

^ ^ 

(pv)  ■ 0 or  PV  - m (50) 

where  m Is  a constant  of  the  motion.  Since  ahead  of  the  wave  v - - 
and  P " Pq  • density  of  the  undisturbed  material, 

pv  - m - - p^Cg  (51) 

Mote  that  this  Integral  of  the  continuous  motion  In  the  steady  wave  takes 
the  same  form  as  the  Jtap  condition  across  a surface  of  discontinuities  In 
the  dependent  variables  (see  [1],  Section  2).  The  reason  for  this  Is  that 
between  any  two  cross-sections  which  propagate  with  the  wave,  an  Invariant 
distribution  of  density,  velocity,  stress  and  stretch  Is  maintained  so  that 
no  contributions  to  the  conservation  relations  are  Introduced  by  the  material 


bstwmen  the  sections 


An  exactly  similar  deduction  applies  for  conservation  of  momentum  and 
energy,  leading  to  the  Integrals  of  the  motion;  Identical  to  the  Jump  rela- 
tions ([1]  Sections  2 and  8) 


-2  „ 
- a + pv  ■ P 


e + V /2  - o v/m  ■ E 


where  e Is  the  specific  Internal  energy,  and  P and  E are  constants. 

In  (53)  adiabatic  loading  Is  considered  as  commonly  acceptable  for  the  ana- 
lysis of  wave  propagation.  Since  no  lateral  strain  occurs  In  plane  waves, 
mass  conservation  demands 


Pq  “ . or  ^ “ Pq/p 


and  (51)  yields 


for  all  X . Thus 


V + XCg  - 0 


[v]  + [X]Cg  - 0 


where  [ ] denotes  the  difference  between  the  values  of  the  enclosed  variable 
at  any  two  cross-sections  propagating  with  the  wave.  Now  since 

[v]  = [v]  (57) 

because  the  two  velocity  variables  differ  only  by  the  superposed  - , 

(56)  becomes,  for  srbltrary  choice  of  axes  in  the  direction  of  wave  propagation: 
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(58) 


Iv]  “ - [A]C^ 

which  is  equlvaleat  to  (38)  when  C • C as  has  been  established  for  the 
case  when  the  material  ahead  of  the  wave  Is  at  rest. 

In  a similar  way  (52)  with  (51)  can  be  written 

-O-  PgCg  V - P 

hence 

(ol  + PoCg[v]  = 0 (59) 

and  using  (58) 

[o]  - PgCg  [X]  (60) 

which  Is  equivalent  to  the  Rayleigh  line  expression  (43).  A similar  mani- 
pulation on  the  energy  balance  equation  (53),  using  the  relation  ([1],  eq. 
(9)) 

[abj  ■ a [b]  + [a]  b (61) 

where  the  superposed  bar  indicates  the  arithmetic  mean  of  the  enclosed 
variable  values,  gives  the  well  known  energy  relation: 

pQlel  - o[X]  - (62) 

Hote  that  if  values  of  the  stresses  and  stretch  ratios  are  known  before 
and  after  passage  of  the  wave,  the  work  done  In  deforming  the  material  and 
the  Increase  In  Internal  energy  Is  given  by  (62)  without  Introducing  addl- 
itional  Information  about  the  constitutive  law  of  the  material. 

In  view  of  the  remarkably  confining  constraints  which  the  existence 
of  a steady  wave  Imposes  on  the  dependent  variables:  stress  o , velocity  v , 
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and  stretch  ratio  X , one  may  wonder  whether  such  waves  are  produced  In 
practice.  That  a close  approximation  to  such  a configuration  Is  achieved 
In  plate  slap  Impact  experiments  has  been  demonstrated  experimentally  by 
Johnson  and  Barker  [11].  Duvall  [12], (p.  102),  gave  a qualitative  descrip- 
tion of  the  stability  of  a steady  wave  In  terms  of  the  distortion  of  the 
Baylelgh  line  and  the  Influence  Introduced  by  this  towards  re-establishing 
the  steady  configuration.  Bland  [13]  Investigated  analytically  the  develop- 
ment of  a steady  wave  for  a Kelvin  type  viscoelastic  solid  with  a nonlinear 
equilibrium  response.  He  showed  by  asymptotic  methods  that  any  monotonlcally 
Increasing  loading  pulse  would  tend  towards  a steady  wave  configuration,  and 
that  for  a step  pulse,  the  time  of  formation  Is  of  the  order  of  five  times  the 
shock  thickness  divided  by  the  excess  convective  velocity.  Particular  cases 
of  elastic-plastic  waves  were  evaluated  by  Clifton  [lA]  who  demonstrated  the 
generation  of  steady  waves.  The  matter  has  been  reviewed  recently  by  Herrmann 
[9]  and  can  clearly  play  an  Important  role  In  the  Investigation  of  rate  effects, 
particularly  at  the  higher  end  of  the  range. 

The  fact  that  the  work  done  on  the  material  as  a shock  wave  traverses  It, 
(62) , Is  not  based  on  the  constitutive  equation  of  the  material  demands  that 
the  equation  adopted,  or  equivalently  the  physical  mechanism  envisaged,  must 
be  flexible  enough  to  permit  the  required  absorption  of  energy.  Part  of  the 
work  Is  expended  In  producing  pure  volume  compression  which  Is  reversible  on 
pressure  release,  the  rest  Is  associated  with  distortion  and  Involves  both 
elastic  and  plastic  deformation.  The  latter  Is  Irreversible  and  Involves 
dissipation  of  mechanical  work  Into  heat  according  to  (23).  A rate-dependent 
or  vlsco-plastlc  law  Incorporates  the  required  flexibility  to  absorb  mechanical 
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energy  since  an  Increase  In  the  rate  of  strain  requires  the  stress  driving 
It  to  Increase  and  hence  also  the  absorption  of  work  for  a given  strain 
Increment.  Since  the  shock-wave  propagation  speed  for  a material  Is  Insensi- 
tive to  changes  In  loading  and  temperature,  for  an  approximately  constant  wave 
amplitude  the  strain  rate  Increases  as  the  wave  becomes  thinner,  and  such  a 
change  In  configuration  can  accommodate  the  energy  absorption  requirement. 

If  the  coefficient  In  the  strain-rate  dependent  expression  for  stress  Is 
large,  only  moderate  thinning  will  be  required  to  provide  a certain  energy 
absorption,  but  a smaller  coefficient  In  the  strain-rate  term  will  call  for  a 
narrower  shock  wave  profile,  and  hence  a higher  strain  rate.  A simple  means 
of  evaluating  this  situation  for  a simpler  elastic-plastic  law  than  that  con- 
sidered In  [1]  has  been  presented  by  Kelly  and  Glllls  [IS]. 

The  rate  Independent  law  suggested  In  [1]  to  permit  the  analysis  of  elas- 
tic-plastic shock  waves  In  a work  hardening  material,  was  to  replace  the  plastic 
work  generation  relation  (20)  by  an  Incremental  relation  across  the  shock: 

pIw]-Y[+eP]  (63) 

where,  as  In  Section  2,  the  square  brackets  Indicate  the  jump  In  the  enclosed 
variable  across  the  shock,  and  the  bar,  the  average  of  the  values  of  the 
m variable  on  the  two  sides  of  the  shock.  It  does  not  appear  that  the  plastic 
work  defined  In  this  way  sould  be  made  consistent  with  the  total  work  defined 
by  (62) , since  there  Is  no  mechanism  to  Increase  the  flow  stress  above  the 
static  yield  to  achieve  the  required  energy  absorption  and  to  meet  the  con- 
straint of  the  Rayleigh  line  relation.  The  fact  that  plastic  flow  must  pro- 
vide the  mechanism  for  the  dissipation  of  energy  was  discussed  In  Section  8 
of  [11. 
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The  steady  wave  analysis  developed  In  that  Section,  based  on  a strain- 
rate  dependent  work-hardening  plasticity  law.  Includes  full  consideration  of 
thermo-mechanical  coupling  according  to  the  elastic-plastic  relations  described 
in  Section  2 of  the  present  report.  Although  Implemented  in  [1]  for  the  hypo- 
thetical rate  law  which  permits  some  plastic  flow  below  the  static  yield,  it 
applies  for  general  rate  laws  which  Incorporate  a strict  static  yield  stress. 

Since  steady  waves  which  propagate  into  material  at  rest  are  normally 
considered,  we  denote  the  wave  speed  based  on  laboratory  or  Euler  coordinates 
by  C for,  as  already  shown, it  is  equal  to  the  Lagrange  wave  speed  defined  on 
the  reference  configuration.  Then,  according  to  (45) , the  stationary  wave 
coordinate  x is  given  by 

X - X - Ct  (64) 


in  terms  of  which  the  stress,  velocity  and  strain  can  be  expressed.  The  materi- 
al derivative  of  a function  g(x)  is  given  by  v(dg/dx) , so  that  the  strain- 
rate  law  (28)  takes  for  form 


V 


(Ej^  - I e^)k(2) 


(65) 


and  the  plastic  work  relation  (29) : 


V ^ ^ (e^  - I eP)\(z)  (66) 

dx  ^0 


The  variable  v can  be  eliminated  from  (65)  and  (66)  by  means  of  the  mass 
conservation  relation  (51),  v ■ m/p  , and  the  conservation  relations  (52) 
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and  (53)  can  be  used  to  provide  expressions  for  the  plastic  strain  and  the 
plastic  work  ({1],  eqs.  (71)  and  (72)): 


(l-Y)w  - E - (e  +^  + 5-  ) - (p_(p+^))(p.(p+^)  + SHE  ) (68) 

where  g(p,0)  Is  the  specific  internal  energy  for  pure  dilatation,  and  is 
associated  with  the  corresponding  free  energy  component  f in  (9).  Eqs. 

(65)  to  (68)  form  a system  of  four  equations  for  the  variables  which  define 
the  state  of  the  material  , 6 , ej  and  w , which  are  all  functions  of 
position  X in  the  stationary  wave.  The  fact  that  two  of  the  equations  are 
differential  and  two  algebraic  gives  a particular  structure  to  the  system 
which  is  discussed  at  some  length  in  Section  8 of  [1].  This  permits  a solu- 
tion in  the  form  of  an  elastic  shock  combined  with  a steady  plastic  wave  in 
which  the  variables  change  continuously.  This  combination  corresponds  to 
the  over-driven  wave  discussed,  for  example,  in  [9],  p.  10.  The  strength  of 
the  wave  can  be  expressed  by  the  mass  flow  m , which  determines  the  elastic 
shock  front  across  which  eP  and  w do  not  change.  The  post  elastic-shock 
values  of  the  variables  then  form  the  initial  conditions  for  solution  of  the 
system  of  equations  which  determines  the  plastic  steady  wave.  Eqs.  (67)  and 
(68)  can  be  used  to  eliminate  eP  and  w from  (65)  and  (66)  which  results 
in  a pair  of  differential  equations  for  and  6 . Integration  from  the  initial 
conditions  is  continued  until  the  rate  of  plastic  strain  drops  to  zero,  and 
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this  condition  determines  the  equilibrium  stress  after  passage  of  the  wave. 

If  the  wave  profile  as  a function  the  space  variable  x Is  not  needed 
but  only  the  ranges  of  stress,  strain  etc.  occurring  In  the  wave,  the  quotient 
of  (66)  and  (65) 


3 P. 

2 


(69) 


can  be  used  to  obtain  a single  differential  equation  which  determines  the 

p 

ranges  of  these  variables.  Substitution  into  (69)  for  and  w from 

(67)  and  (68)  yields  an  ordinary  differential  equation  In  the  (e^^  , 0) 

plane  with  the  post  elastic-shock  values  as  Initial  condition.  The  solution 

values  (c^  , 6)  give  an  Image  of  the  steady  wave  In  this  plane,  and  cal- 

p 

culatlon  of  the  corresponding  values  of  and  w from  (67)  and  (68)  will 

determine  the  extent  of  the  wave  profile  since  It  terminates  when  plastic 
deformation  ceases.  It  will  be  observed  that  In  forming  (69)  as  the  quotient 
of  (66)  and  (65),  the  strain-rate  response  function  (Cj^  “ 3Cj^/2)k(z)  can- 
cels, and  thus  does  not  Influence  the  stress  and  strain  range  of  the  wave. 
This  Is  consistent  with  the  Rayleigh  line  concept  which  Is  not  dependent  on 
the  precise  strain  rate  law,  although  the  space  distribution  of  the  wave  Is. 

For  a smaller  value  of  mass  flow,  m , a precursor  elastic  wave  will  be 
deduced  which  Is  propagated  ahead  of  the  plastic  wave.  The  latter  extends 
from  the  yield  stress,  where  the  plastic  strain  rate  Is  zero  through  a region 
of  plastic  flow  until  an  equilibrium  stress  is  reached  for  which  the  rate  of 
plastic  Flow  Is  again  zero.  Such  a wave  corresponds  to  the  underdriven 
wave  discussed  In  [9],  p.  11. 


These  analyses  incorporated  complete  thermo-mechanical  coupling  according 
to  the  thermo-elastic-plastic  theory  presented  in  Section  2.  The  major  con- 
tributor in  this  regard  is  the  dllatatlonal  component  of  the  internal  energy 
e in  (68),  which  is  perhaps  best  represented  by  the  mle-Gruneisen  equation 
of  state  deduced  from  hydrodynamic  analysis.  Other  analyses  of  the  problem 
have  been  limited  to  the  adoption  of  a purely  mechanical  theory  without  explicit 
thermo-mechanical  coupling.  Clifton's  theory,  [7],  is  limited  to  moderate 
stresses  involving  little  entropy  change,  which  is  assumed  to  be  zero  in  the 
analyses  of  plane  waves  ([17],  pages  135  and  137).  Herrmann  formally  neglects 
thermo-mechanical  coupling  ([9]  p.  2)  but  in  his  analysis  of  plane  waves 
in  an  aluminum  alloy  [16]  he  uses  the  Hugonlot  relation  to  formulate  the 
equation  of  state,  which  Incorporates  thermo-mechanical  coupling  in  the  dila- 
tatlonal  component  of  the  deformation,  which,  as  already  mentioned  normally 
comprises  the  major  component. 

As  mentioned  in  Section  2,  work-hardening  rate-dependent  elastic-plastic 
analysis  introduces  a strain-rate  history  effect,  since  for  a given  strain 
increment,  the  plastic  work  and  hence  hardening  is  greater  at  higher  strain 
rates.  This  is  of  the  nature  of  such  Influences  observed,  for  example,  by 
Klepaczko  [6],  but  may  be  too  pronounced,  since  Herrmann  [16]  cites  evidence 
that  strain  hardening  provides  an  adequate  basis  for  analysis.  However,  the 
structure  of  the  wave  analyses  described  above  including  thermo-mechanical 
coupling  with  an  Internal  variable,  plastic  work,  w , suggests  how  such  a 
variable  can  be  incorporated  into  steady  wave  analysts,  when  it  is  not  tied 
to  the  physical  quantity:  plastic  work.  A functional  relatior  not  so  constrained 
could,  perhaps  be  selected  to  provide  an  adequate  representation  of  the  strain- 
rate  history  influence. 


26 


It  Is  well  known  (Hill  [8]  p.  39)  that  for  rate-independent  elastic- 
plastic  isothermal  theory  based  on  the  Mises  yield  condition,  work  hardening 
and  strain  hardening  give  identical  responses.  For  rate  dependent  analysis 
this  is  not  so,  and  in  particular  the  rate  history  influence  of  work  hard- 
ening already  discussed  does  not  arise  if  strain  hardening  is  assumed.  For 
a prescribed  plastic-strain  Increment  in  monotonlc  loading,  the  equilibrium 
hardening  will  be  Independent  of  the  rate  at  which  the  strain  was  imposed. 

A strain  dependent  version  of  the  theory  already  described  is  therefore  of 
interest.  If  the  complete  thermodynamic  structure  is  maintained  this  change 
will  not  modify  the  theory  appreciably  since  evaluation  of  plastic  work  is 
needed  to  express  the  entropy  Increase  due  to  the  irreversibility  of  plastic 
flow  (23). 

For  a strain  hardening  material,  the  yield  strain  Y in  (18)  would  in 
general  be  a function  of  the  generalized  plastic  strain,  e ^ , and  tempera- 
ture, but  for  monotonlc  loading  in  a plane  wave,  c , so  that: 

Y * Y(ej  , 0)  (70) 

In  the  rate  Independent  theory  applied  to  shock  wave  evaluation,  the  problem 
encountered  with  work  hardening,  that  the  plastic  work  absorbed  as  the  shock 
wave  traverses  an  element  cannot  be  determined  without  analysing  the  shock 
structure  still  arises,  for  plastic  work  appears  in  the  energy  conservation 
equation.  A rate  of  strain  term  is  needed  to  determine  the  shock  structure. 
Similarly  in  the  plastic  steady  wave  structure  evaluation,  equations  (65)  to 
(68)  are  still  all  needed,  since  plastic  work  appears  on  the  left  hand  side 
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of  (68).  Thus  Che  fora  of  solution  Is  the  same  for  a strain-hardening  material 
as  for  a work-hardening  one. 

Ho%rever,  the  term  (1  - y)v  on  the  left  hand  side  of  the  energy  balance 
equation  (68)  is  likely  to  be  small  compared  with  other  terms  in  the  equation 
if  Che  impact  stress  is  many  times  larger  than  the  yield  stress  which  is  usu- 
ally the  case  in  place  slap  experiments.  The  term  (1  - y)v  represents  the 
part  of  Che  plastic  work  which  is  scored  in  Che  material  as  dislocations  and 
lattice  ucfects  following  plastic  flow.  The  factor  (l  - y)  'v  0.1  , for  most 
of  the  plastic  work  is  dissipated  directly  into  heat.  The  plastic  work  it- 
self will  be  much  smaller  chan  other  terms  in  the  energy  balance  equation 
since,  for  example,  it  is  given  by  the  yield  stress  times  the  plastic  strain 
increment  while  Che  energy  absorbed  in  volume  compression  is  of  the  order  of 
the  total  stress  multiplld  by  the  strain  increment.  Thus  very  little  error 
will  be  Involved  in  neglecting  the  term  (1  - y)v  and  then  Che  analysis 
simplifies  since  plastic  work  does  not  appear  as  a variable.  Eq.  (68)  then 
gives  Che  image  of  the  shock  structure  in  the  (e^,e)  plane  and  (67)  determines 
the  plastic  strain.  The  differential  relation  (65)  is  needed  only  to  obtain 
the  wave  profile  in  space.  Thus  a much  simpler  solution  for  the  steady  wave 
results.  In  view  of  this  simple  structure,  for  practical  calculations  it 
would  be  advantageous  to  consider  waves  of  increasing  and  decreasing  strain 
separately  and  use  a direct  overstress  form  for  the  plastic  strain  rate  in- 
stead of  (28)  in  terms  of  Che  vsrlable  z , which,  because  of  its  quadratic 

p 

nature  applies  to  both  type  of  waves.  In  that  case  the  factor  (c^  - 3e^  / 2) 

•P 

would  not  be  needed  to  generate  the  appropriate  sign  for  , which  could  be 
included  separately. 


4. 


Discussion  and  Conclusions 


The  analysis  of  shock  wave  propagation  through  elastic-plastic  materials 
which  exhibit  work-hardening,  strain-rate  dependent  response  to  stress  requires 
the  solution  of  a differential  equation  associated  with  the  structure  of  a 
steady  plastic  wave.  The  theory  Includes  thermo-mechanical  coupling,  and  a- 
bout  lOZ  of  the  plastic  work  Is  assumed  absorbed  in  the  resulting  dislocation 
and  lattice-defect  distribution  as  measured  by  Qulnney  and  Taylor  [17].  This 
component  of  the  plastic  work  absorbed  Is  small  compared  with  the  total  work 
done  as  the  wave  traverses  a section,  so  that  a negligible  effect  on  the 
process  will  arise  from  assuming  this  component  of  the  absorbed  energy  to  be 
dissipated  Into  heat  as  Is  the  rest  of  the  plastic  work.  In  the  case  of  strain- 
hardening, this  minor  change  in  the  thermo-mechanical  coupling  process  permits 
the  stress  generated  across  a plastic  shock  wave  to  be  evaluated  using  the 
conservation  relations  only,  which  are  algebraic,  rather  than  differential  In 
nature.  This  offers  the  possibility  of  much  more  convenient  evaluation  of  the 
stress  range  and  total  and  plastic  strain  generated.  For  rate  dependent 
elastic-plastic  materials  work-hardening  and  strain-hardening  laws  lead  to 
qualitatively  different  behavior,  which  Is  not  so  for  the  corresponding  rste- 
Independent  theory.  Strain-hardening  ^^ppoars  to  reproduce  exoerlmental  flndlnss 
and  has  been  used  by  Herrmann  In  his  studies  of  wave  propagation  [9,  16]. 

Since  the  relations  are  algebraic,  the  Jumps  In  the  dependent  variables  can 
be  determined  without  evaluating  the  shock  structure,  and  it  may  be  adequate 
In  many  cases  to  utilize  this  Information  In  design  calculations. 

Much  of  the  experimental  work  on  the  measurement  of  stress-wave  profiles 
has  been  concerned  with  elucidating  the  strain-rate  response  law  governing 


plastic  flow.  As  described  In  Section  3,  the  relationship  between  wave 
profiles  and  the  rate  of  strain  law  Is  associated  with  the  Increase  In  stress 
needed  to  lift  the  equilibrium  stress-strain  relation  up  to  the  dynamically 
Imposed  Rayleigh  line  Including  the  thermal  Influences  on  these  relations. 

However,  possible  difficulties  of  Interpretation  arise  In  the  case  of  materials 
having  a marked  microstructure  since  this  Itself  will  also  tend  to  broaden 
the  wave  front  and  thus  generate  a lower  average  strain-rate.  In  the  case  of 
poly-crystalllne  metals,  this  can  arise  from  anisotropy  of  the  crystallites, 
and  particularly  in  the  case  of  Initiation  of  plastic  flow,  because  of  the 
concomitant  residual  stress  In  the  crystallites  resulting  from  the  process 
used  to  form  the  specimen.  These  Influences  have  been  commented  on  by  Herrmann 
[9].  For  porous  media  the  effect  of  Inhomogenelty  can  be  much  more  marked, 
and  the  local  micro-motion  of  the  collapsing  pores  for  an  element  In  com- 
pression can  have  an  appreciable  influence  on  the  average  stress-strain  rela- 
tion. It  was  shown  in  [18]  that  the  inertia  forces  of  the  material  adjacent 
to  the  collapsing  pores  necessitate  an  Increase  In  the  average  compressive  | 

stress  sccordlng  tc  a strain- rate  dependent  relation.  This  Is  so  since  | 

even  at  constant  average  rate  of  strain,  the  surface  of  the  closing  pores  must 
accelerate  as  the  pore  surface  decreases  In  area  In  order  to  accoimodate  the 
effectively  incompressible  deformation  due  to  plastic  flow.  This  requires 
additional  compressive  stress  proportional  to  the  strain  rate.  This  con- 
stitutes a system  rate  of  strain  Influence  superposed  on  whatever  material  I 

rate  of  strain  characteristic  exists.  Such  an  effect  will  contribute  to  the  j 

wave  profile  generated.  It  may  therefore  be  Important  to  Incorporate  such  | 

effects  Into  investigations  of  basic  material  strain-rate  Influences.  Moreover  j 
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the  existence  of  Inhomogeneity  will  modify  wave  profiles  evaluated  on  the 
assumption  of  homogeneous  material,  which  may  dictate  against  placing  reliance 
on  computed  wave  profiles.  This  may  lend  additional  encouragement  to  attempts 
to  satisfy  design  requirements  on  the  basis  of  calculations  of  jumps  of  stress 
and  strain  across  the  wave,  while  reducing  the  significance  of  the  detailed 
stress  profile  to  a minimum. 
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